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We present a first calculation of transverse momentum dependent nucleon observables in dynam- 
ical lattice QCD employing non-local operators with staple-shaped, "process-dependent" Wilson 
lines. The use of staple-shaped Wilson lines allows us to link lattice simulations to TMD effects de- 
termined from experiment, and in particular to access non-universal, naively time-reversal odd TMD 
observables. We present and discuss results for the generalized Sivers and Boer-Mulders transverse 
momentum shifts for the SIDIS and DY cases. The effect of staple-shaped Wilson lines on T-even 
observables is studied for the generalized tensor charge and a generalized transverse shift related 
to the worm gear function giT- We emphasize the dependence of these observables on the staple 
extent and the CoUins-Soper evolution parameter. Our numerical calculations use an n/ = 2-1-1 
mixed action scheme with domain wall valence fermions on an Asqtad sea and pion masses 369 MeV 
as well as 518 MeV. 
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The picture of the nucleon as a system of interacting quarks and giuons naturally leads to the question about 
the intrinsic motion of these elementary particles inside the proton or neutron. This intrinsic motion, specifically 
with respect to the transverse momentum, can be described in terms of Transverse Momentum Dependent Parton 
Distribution Functions (TMDs), see, e.g., chapter 2 of Ref. [T] for a recent review. TMDs for quarks, generically 
denoted by fi{x,k'^), gi{x,k'^), etc., encode essential information about the distribution of partons with respect to 
the longitudinal momentum fraction, x, and intrinsic quark transverse momentum, fcx- With certain restrictions in 
mind, they have an intuitively appealing interpretation as three-dimensional probability densities (21 [3]. TMDs can, 
for example, be studied on the basis of angular asymmetries observed in processes such as Semi Inclusive Deep Inelastic 
Scattering (SIDIS) using suitable QCD factorization theorems that go beyond the standard collinear factorization, see, 
e.g., Refs. [IHl]. In contrast to the usual collinear PDFs, TMDs turn out to be in general non-universal, i.e., process- 
dependent. The process dependence arises from the difference in the final and initial state interactions in SIDIS 
and Drell-Yan scattering, respectively. On the theoretical level, it can be understood as an intriguing consequence 
of the local color gauge invariance of the strong interaction and the corresponding non-trivial gauge-link structures. 
Specifically, QCD factorization leads to the remarkable prediction that the naively time reversal odd (T-odd) TMDs, in 
particular the Sivers and Boer-Mulders functions, differ in sign for DY compared to SIDIS, jT-odd,siDis _ _yT-odd,DY^ 
The implications and consequences of these observations continue to stir intense interest of many theoreticians and 
experimentalists, as a number of fundamental questions and interesting puzzles remain to be addressed. Motivated 
by promising experimental results from COMPASS, HERMES and JLab (see, e.g., [SHTU] and references therein), as 
well as considerable progress on the theoretical and phenomenological sides during recent years, an essential part of 
the physics program of future facilities will therefore be targeted in this direction, including JLab 12 GeV and the 
proposed EIC at JLab or BNL. 

Theoretical calculations of TMDs from first principles require non-perturbative methods such as lattice QCD. In 
previous works, we have introduced and explored techniques that allow the computation of the underlying amplitudes 
on the lattice using non-local operators [TTlfTS] . Our numerical studies for a "process- independent" , direct gauge link 
geometry already produced encouraging results. In this work, we present a first exploratory lattice study employing 
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a more complex, "process-dependent" link geometry that gives us rather direct access to highly interesting T-odd 
observables. 

In section [nj we present the formalism and techniques required for our calculations, and provide definitions of the 
relevant T-odd and T-even TMD observables. After a short introduction to the lattice computations at the beginning 
of section |III[ we continue with a presentation and discussion of our numerical results for the generalized shifts and 
tensor charge. A summary and conclusions are given in section [TV] 



II. FORMALISM 



A. Definition of TMDs 



In a relativistic quantum field theory, the question "What is the probability to find a quark with a given momentum 
k inside the proton?" needs to be stated more precisely. First of all, it turns out to be advantageous to formulate 
everything in light cone coordinates, see appendix [Aj and to consider a frame of reference where the nucleon has large 
momentum in z-direction, i.e., P+ ^ mj^, Pt — 0. In light cone coordinates, the components fc"*", fcx, k~ of the 
quark momentum k scale as /rriN, 1, ttin/P^, respectively, under boosts along the z-axis. Thus the longitudinal 
momentum fraction of the quark x — k^ / P'^ and its transverse momentum fcx are invariant under boosts along the 
z-axis, while the k" component is suppressed. This leads to the concept of transverse momentum dependent parton 
distribution functions (TMDs), which are functions of the longitudinal momentum fraction x = k'^ jP^ and of the 
quark transverse momentum fcx- The transverse momentum components fex are particularly interesting, because 
they describe an intrinsic motion of the quarks inside the proton that occurs independent of the momentum of the 
proton itself. This gives us a unique picture of the dynamics inside the proton. Moreover, the TMDs are an important 
ingredient in our understanding of the origin of large angular- and spin-asymmetries found in experiments studying, 
e.g., semi-inclusive deep inelastic scattering (SIDIS) or the Drell-Yan process (DY). 

In a naive approach based on a theory quantized on the light front, one obtains a momentum dependent number 
density of quarks from fcx) ^ \ X]a=±i Sa=±i \'^\,qk^i ^t) \P, S) p (up to normalization factors), where a\^q is 
an annihilation operator of quarks of flavor q and helicity A. The average over nucleon helicities | X]a±i implements 
an average over the spin S in the nucleon state \P, S). In this example, the TMD /i(x, fcx) describes the distribution 
of unpolarized quarks in an unpolarized nucleon. Rewriting the annihilation operator in terms of local quark field 
operators q and q reveals a problem: fi{x, fcx) is a Fourier transform of the matrix element (P, S\ q(0)^~^ q{b) |P, S) with 
respect to the position b, and the bi-local operator g(0)7+g(5) is not gauge invariant, see Ref. [J for a review of the 
issue. Gauge invariance can be restored by inserting a Wilson line U [Cf,] between the quark fields, as defined in appendix 
[A| The Wilson line introduces divergences that cannot be treated by conventional dimensional regularization |14j . 
Several different schemes have been proposed in the literature as to how to subtract those divergences [3l El |6l [7l [14]- 
flQ] . see Ref. [20^ for a recent comparison. In general, these schemes require the introduction of a so-called soft factor 
iS inside the defining correlator of TMDs. The starting point for our discussion of TMDs is thus a correlator of the 
general form 



$[i^l(fc,P,5;...) = /T^e^'^''^ ^ (1) 

^ ^ J (27r)4 5(&2;...) ^ ' 

The detailed properties of the Wilson line U[Ch] and the soft factor need to be specified by additional parameters, 
which we indicate by the dots ". . ." for now and which will be discussed later. Moreover, all objects above implicitly 
depend on a UV renormalization scale /i. 

In Eq. ([T]) , S stands somewhat symbolically for an expression that can, depending on the formalism, involve several 
vacuum expectation values. For example, in the scheme developed in Refs. O III [19], our factor 5(6^; . . .) would be 
(using the notation of those references) 



S{b';. 



where each of the objects 5(0) (^t, • ■ ■) is a vacuum expectation value of Wilson line structures. In this specific 
framework, the starting point of the discussion is space-like Wilson lines. Some Wilson lines remain tilted away from 
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the light cone, leading to the dependence on the rapidity parameter in the above expression. Other Wilson lines, 
including those contained in U[Ch\ in the numerator of Eq. ([T]), are brought back to the light cone in the sense of a 
limit, as indicated symbolically by +00 and —00 in the equation above, see Refs. [SIE] for details. As we will see in 
section [II E| and below, certain matrix elements with space-like structures of Wilson lines are directly accessible on the 
Euclidean lattice, whereas taking the light-cone limit is only possible in the form of a numeric limit and technically 
challenging. For the purposes of our treatment, however, we do not need to go into any detail concerning the definition 
of 5(5^; . . .) in any particular framework, since it will cancel in the observables we consider. 
Integrating the correlator over the suppressed momentum component k~ yields 



$1 



d{b-P 



ix{b-P)-ibT-kT 2^ 



\ (P,5| g(0)rz^[a] q{h) |P,5) 



(27r)2 J {2n)P+ 



(3) 



b+=0 



Notice that integrating over corresponds to setting 6+ = 0. As a consequence, x (b-P) and fex <-> bx act as 
independent pairs of Fourier conjugate variables in the expression above. The above correlator can be decomposed 
into TMDs. For choices of the Dirac matrix F that project onto leading twist, one obtains 



<i>[^+l(a;,/cT;P,^,...) 
$[^^^'l(a;,fcT;P,5,...) 



gjj kj Sj _|_ 
/it 



= A51 H giT , 



= hi 



rriN 

S^r "it 



(4) 
(5) 
(6) 



The TMDs /i, (71, hi, giT-, hiL^ ^iti fir ^^'^ ^1 functions of x, fc^, ^ and further parameters related to 
regularization and link geometry. The structures shown in brackets [ ]odd involve so-called naively time reversal odd 
(T-odd) TMDs, namely the Sivers function [2F and the Boer-Mulders function h-^ [55]. The origin of the above 
paramctrization and the special role of T-odd TMDs will become clear after we have discussed the geometry of the 
gauge link path Cf, and symmetry transformation properties. 



B. General strategy 

At this point, several remarks are in order as to how we aim to introduce TMD observables that can be accessed 
with lattice QCD using a non-local operator technique. Due to the underlying operator structure, the situation is 
quite different from that of standard coUinear PDFs and offers unique opportunities and challenges. 

As an introductory example, consider the definition of a standard PDF in the unpolarized case, 

/l(x)= ^ / ^^"e"^^'" (^'^1 mi^m,nh-] q{nb-) \P,S) . 

For PDFs, the gauge lmkU[0, nb^] is simply a straight, light-like Wilson line of finite extent connecting the two quark 
field operators [27]. No continuous Lorentz transformation exists that allows us to "rotate" the non-local operator 
q{0) 7+ U[0, nb~] q{nb~) into Euclidean space. The light-like separation stays always light-like, but in Euclidean space 
objects cannot have any extent in (Minkowski-) time. As a consequence, one is forced to invoke the operator product 
expansion to cast the calculation in terms of local matrix elements which can be accessed using lattice QCD. 
The situation for TMDs differs fundamentally in several aspects: 

1. The separation b of the quark field operators has an additional transverse component, b = nb^ + b±. Thus, 
in general, this separation is space-like. This opens the possibility of a direct representation of the non-local 
operator in Euclidean space. 

2. The geometry of the gauge link U[Cb] is more complicated, depends to a certain degree on the experiment under 
consideration and in general extends out to infinity. As a result, it becomes questionable whether an expansion 
in terms of local operators is possible at all. 
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3. Regularization is more complicated, leading to the introduction of the soft factor S and additional regularization 
parameters beyond the usual renormalization scale fi of the MS scheme. 

The first two items listed above are our main motivation to develop a technique for lattice studies of TMDs based 
on non-local operators. It should be emphasized that this technique can only work for the analysis of certain TMD- 
related observables within a limited kinematical range. The method cannot be applied to study the x-dependence of 
PDFs directly, without the use of non-trivial extrapolations. 

In previous publications [12|, I13j . it was demonstrated that the non-local operator technique is quite promising 
and produces interesting results, for a simplified gauge link geometry, at least on a qualitative level. The crucial 
connection between the formalism in Minkowski space and the results from Euclidean space is provided through a 
parametrization in terms of invariant amplitudes^ Ai{b'^ ,b-P). By virtue of their Lorentz-invariance, the calculation 
of these amplitudes can be performed in any desired Lorentz frame. In particular, for the generic off-light cone 
kinematics appropriate for TMDs, there is no obstacle to performing the calculation in a frame in which the nonlocal 
operator in question is defined entirely at one fixed time. In this frame, one can cast the computation of the nonlocal 
matrix element in terms of a Euclidean path integral, evaluated employing the standard methods of lattice QCD. 

The study at hand builds directly on Ref. |13| . and we refer the reader to that publication for an introduction 
to the essential principles of the methology. One of the remaining challenges identified in Ref. [13] concerns the 
geometry of the gauge link. In the present study, we replace the simple straight connection by a staple-like path that 
corresponds more accurately to the situation in phenomenology. We stress that these gauge link structures are part 
of the established phenomenological framework, which we take as given, and not a new assumption related to our use 
of lattice QCD as a calculational method. Whereas our results depend on the gauge hnk structure, specific physical 
processes such as SIDIS and DY unambiguously correspond to definite instances of that structure. Throughout our 
discussion, we clearly identify the SIDIS and DY limits of our data. 

It is important to point out that our assumptions about the operator structure of TMDs rely on factorization 
arguments that are much more involved than for the usual PDFs. In fact, one must be judicious concerning the classes 
of reactions for which it can be assumed that a factorization framework with well-defined TMDs exists. For example, 
it has been realized recently [28, 29J that TMD factorization generally fails for large reaction classes, in particular 
processes with multiple hadrons in both the initial and the final state. While the consequences of this observation 
are not yet all known, it appears certain that to develop a TMD framework for these processes, a fundamental 
change of perturbative QCD techniques is needed. It could, e.g., very well be that measurement-independent cross- 
sections are simply not defined for certain reaction classes and that, instead, the appropriate quantities will be 
entanglement amplitudes which then have to be folded with quantities encoding the measurement process |30j . In 
contradistinction to the aforementioned classes of reactions, for other types of processes such as SIDIS and DY, recent 
progress [6j [71 [19] indicates that a valid definition of TMDs based on factorization arguments indeed is possible, 
within a scheme regularized employing space-like links. Promising steps have been taken to develop the predictive 
capabilities of this framework |31j . A pertinent discussion is given in section I of our previous publication [!13 , with 
further details to be found in the references therein and recent overviews in Refs. [201 [H]- The point which we wish to 
emphasize here is that it is not the purpose of our present work to critique or justify the various approaches to defining 
TMDs in terms of operators and matrix elements which have been advanced in response to issues of factorization and 
regularization. Instead, we will assume that a good definition of TMDs with a connection to phenomenology through 
a valid factorization argument exists for certain classes of processes such as SIDIS and DY; we focus exclusively on 
those TMDs and do not aim to contribute to discussions of factorization, fragmentation functions, or related matters. 
Our starting point thus is the definition of TMDs in terms of a TMD correlator of the rather general form ([T]). 

Working from this definition, we will moreover restrict ourselves to observables in which the soft factor cancels, so 
that specifics of the soft factor are not relevant for our results. Regarding the necessary regularization of the gauge 
link, we pick the proposal that is most suitable for our purposes, namely tilting the gauge link slightly away from the 
light cone i32J, in a space-like direction [HIIZI. We stress that, in choosing this approach to defining and regularizing 
TMDs, we are led to consider kinematics off the light cone from the very beginning, which makes a connection to 
Euclidean lattice QCD feasible, as already noted further above. We emphasize that, within this work, we do not aim 
to arrive at any statements concerning the formal nature of the light-cone limit. We will, however, focus particularly 
on the behavior of our numerical results as we extend the kinematic region as far towards the light cone as possible. 

One necessary step of a lattice calculation is to discretize the operators. The discretization of non-local operators 
as we encounter them here is still a rather new concept. An important assumption we make is that non-local lattice 
operators composed of structures much larger than the lattice spacing essentially renormalize in the same fashion as 



Note that the symbol I in Ref. |12l I13| corresponds to —6 in the present study. 
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their counterparts in the continuum, except that the renormahzation parameters are specific to the lattice action and 
the discretization prescription. We have given reasons for this assumption and explored it numerically in Ref. jl3j . 
see in particular sections III D, IV B, IV C and appendices B, D, G and H therein. However, we point out that a 
more rigorous treatment would still be desirable. Especially the question of mixing properties as one attempts to 
make contact with the local operator formalism remains a challenge for the future. 

Keeping the above remarks in mind, it is worthwhile summarizing the logic underlying our treatment succinctly 
before laying out the details further below: 

1. We start from a definition of TMDs in terms of the correlator ([T]), considering generic off- light cone kinematics 
from the very beginning. 



2. The correlator ([ij is parametrized in terms of Lorentz- invariant amplitudes, cf. section II F below. This crucial 
step permits one to transform results into different Lorentz frames in a simple manner. 

3. On this basis, we choose the Lorentz frame in which the nonlocal operator entering ([T]) is defined at one single 
time as the one most suitable for our calculation. We stress again that there is no obstacle to this choice, since 
the separations in the operator are all space-like. 

4. In the aforementioned frame, the computation of the nonlocal matrix element can be cast in terms of a Euclidean 
path integral and performed employing the standard methods of lattice QCD. 

5. We form appropriate ratios of the extracted invariant amplitudes in which soft factors and multiplicative renor- 
mahzation factors cancel, such as the fcT-shifts discussed in section [lI D[ which, in principle, represent measurable 
quantities. We particularly study the approach to the SIDIS and DY limits in these quantities. 



C. TMDs in Fourier space and a;- integration 

In essence, the lattice method we use allows us to evaluate the ^-dependent matrix elements ^I^g^i^t^ {b,P,S; . . .) 
introduced in Eq. Q. As a result, it is more direct and natural to state our results in terms of Fourier-transformed, 
bx-dependent TMDs and their bx-derivatives. For a generic TMD / we define 

;...) = I d^kre'"^"^ f{x,k^^;...)^27r J d\kT\\kT\ M\bT\\kT\) f , k^T, ■ ■ ■) , (7) 
/(")(x,b^..)^n!f-^9,.y />,b2;...) = ^V^ /d|fcT||feT|f|^yj„(|bT||fcT|)/(x,4;---), (8) 

where the J„ are Bessel functions of the first kind, and mAr is the mass of the target hadron. These objects and their 
potential phenomenological relevance have been discussed in detail in Ref. |33) . Moreover, evolution equations are 
naturally expressed in terms of the Z*^"-*, compare, e.g., Ref. [34] . In the limit [bxl — >■ 0, one recovers conventional 
fcx-moments of TMDs: 

/(")(x,0;...)=|d2fcT(^2^) /(x,fe2;...)EE/(")(.T) . (9) 

However, it is known |35| that fcx-moments like ff^ {x) and /j^^^"* {x) are ill-defined without further regularization. 
The problem is that the integral in the above equation diverges if the integrand does not fall off quickly enough in 
the region of large fcx, where the TMDs /(x, fc^) are perturbatively predictable. Even though fcx-nioments may be 
more familiar to the reader, we therefore do not attempt to extrapolate to bx — 0, but rather state our results at 
finite |bx|, where the fcx-integrals of Eqs. Q and ([8| can be shown to be convergent in the relevant cases [33] . 

Information about the x-dependence of TMDs can be obtained from the lattice via the Fourier-conjugate variable, 
b-P jlHITni. However, the calculations performed in Euclidean space only allow us to access a limited range of 6-P, 
precluding us from performing a straightforward Fourier transform. In this work, we limit ourselves to the study of 
x-integrated TMDs 

;...) = y\x/(x,fe^;...) . (10) 

These are accessible from the data at h-P ~ 0. Here, the superscript denotes the first McUin moment in x. The 
integration is performed over the full range of x. TMDs evaluated at negative values of x can be related to anti-quark 
distributions, see, e.g., [T31[23] for details. 
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D. Quantities suitable for lattice extraction 



Certain ratios of fcT-moments of TMDs have interesting physical interpretations. For example, consider 



rriN 



(11) 
(1,0) 



where 7+ projects on leading- twist. In the context of the density interpretation of TMDs mentioned in section II A 
the ratio above yields the average transverse momentum in ^/-direction, for quarks with given longitudinal momentum 
fraction x inside a proton polarized in x-direction. We will show below that quantities like this can be calculated 
rather directly on the lattice. For the reasons mentioned above, we limit ourselves to ratios formed from x-integrated 
quantities. Let us therefore consider 

j-L[l](l) 

{ky)Tu = mN^^^ . (12) 
Jl 

Ignoring the role of anti-quarks, this ratio, called in the following "Sivers shift" , represents the average transverse 
momentum of unpolarized ("U") quarks orthogonal to the transverse ("T") spin of the nucleon. Note, however, that 

the denominator /j^'^*^^ arises from a difference of quarks and anti-quarks and thus gives the number of valence quarks 
in the nucleon. On the other hand, in the numerator f^^^^^\ the average transverse momentum of quarks and anti- 
quarks is summed over |13[ I23j. A profound interpretation of /j^"'^'*'"'^^ in impact parameter space has been given in 
Ref. ^36j|- However, as mentioned before, understanding /^y^'^^'* simply as a fcx-weighted TMD is problematic, since 
the fcx-integral is expected to be UV divergent. A natural way of circumventing this divergence is to generalize the 
Sivers shift to an expression in terms of the Fourier-transformed TMDs: 



fJ-[l](l)^.2 . N 

(fc,)™(b^; . . .) ^ Zm^2^'''' ■ (13) 

Jl ("T, • ■ •) 



This is the type of quantity that we investigate in the present study. In the limit 6^ = we recover the Sivers 



shift (12), because the Fourier transformed TMDs f^r^ and /} coincide with the moments f^r^ and /{ 



respectively. We are, however, interested in the generalized Sivers shift for non-zero b^, where the said UV-divergence 
disappears. The variable effectively acts as a regulator. Moreover, the bx-dependence allows us to study differences 
in the widths of distributions on a qualitative level. 



E. Link geometry 

The prescription for the geometry of the gauge link path Cf, affects both the number of allowed structures appearing 
in Eqs. (|4])-([6]) and the numerical result for the TMDs. We therefore need to ask which link geometries are appropriate. 

The simplest link geometry is a straight line connecting the quark fields at and b, see Fig. [Ta[ TMDs with 
straight gauge links have been studied on the lattice in Refs. [TH|T3]. While these "process- independent" TMDs are 
interesting from a theoretical point of view in their own right, it is so far not known how to relate these quantitatively 
to the TMDs that play a role in scattering experiments. The operator with straight gauge links offers the largest 
possible degree of symmetry. As a result, T-odd TMDs vanish for straight gauge links. 

For TMDs that allow us to describe measurable effects in scattering experiments such as SIDIS or DY, the form of 
the gauge link is largely dictated by the physical process. To understand scattering experiments at high momentum 
transfer Q, one tries to apply approximations valid for large Q that separate hard, perturbative and soft, non- 
perturbative scales in the dominant physical processes in order to arrive at an expression for the cross section in 
factorized form. In the standard coUinear approximation, all internal transverse momenta are integrated out and 
conventional parton distribution functions and fragmentation functions are used to describe the process. In certain 
kinematical regions this approximation is insufficient. An example is SIDIS, where the momentum P/j of one of 
the final state hadrons is measured after a lepton- nucleon collision at large momentum transfer Q. The transverse 
momentum dependent formalism is needed when the transverse momentum component Phi_ is small with respect to 
Q, see, e.g., Ref. |35j for an in-depth discussion. 

The leading diagram for SIDIS is shown in a simplified, factorized form in Fig. [2] The lower shaded bubble in 
the diagram represents the structure parametrized by TMDs. A gauge link in the TMD correlator arises naturally as 
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an idealized, effective, resummed description of tlie gluon exchanges between tlie ejected quark and tlie remainder of 
the nucleon in the evolving final state, see, e.g., Ref. [37^ for a review. The gauge link roughly follows the direction 
of the ejected quark, in SIDIS by convention denoted by the light-cone n direction. The TMD correlator obtained 
from the squared amplitude thus has parallel Wilson lines attached to each of the quark field operators at and b, 



extending out to infinity along a direction v k, n, see Fig. lb Due to the fact that the gauge link is only an effective 



representation of final state interactions within a framework of suitable approximations, there is a certain degree of 
freedom with respect to its geometry, in particular with regard to the choice of its direction v. At tree level, the most 
convenient choice is an exactly light-like gauge link, v = n. However, going beyond tree-level, it has been found that 
the light-like link introduces so-called rapidity divergences that are hard to remove, see Ref. [38] for a review. One 
way of regulating these divergences is to use a gauge link slightly off the light cone [17], see Refs. [HOI?] for the 
application to SIDIS. In Ref. [3], the direction v is chosen time- like. More recent work in Refs. [51 [7] is based on 
space-like Wilson lines, motivated by the insight that TMDs with this choice of link directions feature a "modified 
universality", i.e., they are predicted to be numerically equal for both SIDIS and DY [T5| up to the expected sign 
changes of T-odd TMDs. The space-like choice of Wilson lines also opens up the possibility of implementing the gauge 
link directly in lattice QCD. 

In Fig. |lb[ the two parallel Wilson lines are connected at the far end by another straight Wilson line. The complete 
gauge link thus has a staple-like shape. Bridging the transverse gap is necessary to render the operator gauge invariant 
and proves to be essential if the light-cone gauge n ■ A — is used [321 SO]- In a covariant gauge, the connecting link 
at infinity can be omitted; this has been exploited in Refs. [HIS]- Lattice calculations are typically performed without 
any gauge fixing. We therefore prefer the notation with an explicitly gauge invariant operator. Moreover, in our study 
we take the limit of an infinite "staple extent" rj explicitly. The gauge link employed in this work thus reads 

ll[ci'"''>]^U[0,vv,Tjv + b,b] , (14) 

where v is space-like. Even at finite rj, this gauge link geometry fulfills the desired symmetry transformation rules, 
as listed in Eq. (C6) of Ref. [13] and discussed further below. Here, we will be mostly concerned with the lowest 
a;-moment of TMDs, corresponding to the case b~ = b'^ = 0. In this case, the connection at the far end is purely 
transverse. 
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We choose v space-like and, as in Refs. [HO [7], we consider TMDs for the choice vt = 0. The Lorentz-invariant 
quantity characterizing the direction of v is the parameter = 2v-P/ y/\v'^\. The hght-hke direction v = n can be 
approached in the limit C ~^ oo- The parameter ^ can be understood as an artificial scale or cutoff introduced to 
regulate rapidity divergences. Within their work on e+e~ -scattering, Collins and Soper provided evolution equations 
for the dependence on ( applicable for C ^ J^qcd [22] ■ Similar equations have been worked out for all leading- twist 
and spin-dependent parton distributions [31j based on the formalism of Ref . [J] . For the more recent formalism of 
Refs. [H[7], evolution equations are presently available for the unpolarized case and the Sivers function [IH]. The 
vectors P and v can be written in terms of rapidities yp and respectively: = vnj^e^^^ j \f2 and v'^ /v~ = — e^^". 
Rewriting ^ as a dimensionless quantity, 

V ■ P 

C = C/2mAr = — =— = ==sinh(j/p-y„), (15) 

I V -T 

reveals that it is essentially a rapidity difference. Notice that the entire system can always be boosted to a frame 
where v has only spatial components, w*^ = 0, = 0. This is crucial for the lattice approach. 



F. Parametrization of the correlator 



The translation of our results obtained in Euclidean space into TMDs defined and interpreted in the con- 
text of light cone coordinates is mediated through a parametrization of the correlator $unsubtr. in terms of 
manifestly Lorentz-invariant amplitudes. For our purposes, it will be important to take the dependence on 
the link direction v explicitly into account. A parametrization of the correlator $unsubtr.(fc, -P, "S*; oou, /i) = 

/ d'^bl{2'iT)'^e^^-^\ {P,S\ g(0)r Ulc'^^'^qib) \P,S) for link paths that extend to infinity into a direction v has 
been worked out in Ref. [23] and involves 32 independent amplitudes Ai and Bi that depend on the Lorentz-invariant 
quantities fc^, k ■ P, k-v/v-P and C. Appendix C of Ref. [33] shows that a parametrization of the corresponding 
6-dependent correlator 'I'unsubtr. is of the same form as the parametrization of $unsubtr. if we substitute k — > —im?j^b. 
We thus obtain 



ubt. - ^nM - '^e^'^P^P.b^v.S^B, (16) 
^^!7nlbt. - <{b-S)A, + '-^{vS)B, (17) 
^^I7nlbtr. ^P^A,^ imlb^ I3 - im^e^^'^'^^P^b^Sp A,, 

2 - 3 



{vpy ^ ' v-p 



- ^(6.5)e^''"^P.6„z;0 Pg - —^{v-S)e^^'^P P^b^vpB^^ (18) 

v-p (v-py 



^^unsubtr. - -mAT^Me + imN{b-S)P^' Aj + m%{b-S)b^As 
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^ ib-Sy'B,3 ~ j^ivSy'^B.i (19) 



vP ^ ^' " vp-- ' y.p 

— —[b-b)v' Bi3 - - — 

V-P [V-P)-' 



I'^tnZbtl = iniNe^^'^P^Ppb^Ai + pl^^-'Ug - im^fcl^^'^Uio - m%{b-S)PH''^ An 

- '^e'^'^'-'PpV^B^ + '-^e^-^P-b.v.B, + '^v^^S'^^B,, ~ '-II^ib-S)PK'^] B,, 

V-P V-P V-P V-P 

- '^{b-S)bK-^^B,r - '-^{v-S)pH-^^B,s + ^{v-S)pK-^^B,, - ^iv-S)bK-^^B,o , (20) 

V-P V-P V-P (V-P)^ 

where b^ . The structures above are compatible with the transformation properties of the correlator 

under the symmetries of QCD. For completeness we list them again in appendix [B] 

Our previous studies of TMDs on the lattice [TU [TJ] were carried out with straight gauge links. In that case only 
the T-even structures involving amplitudes of type Ai appear in the parametrization. As pointed out already in 
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those references, there is not necessarily a one-to-one correspondence between the Ai and Ai (or the Bi and Bi). For 
example, As contributes to Aq, Aj and Aq. Note that I = —b in Refs. [T^ [T5] . 

In the above parametrization, factors of {v-P)~" ensure that the structures are invariant under rescaling of v, i.e., 
V — )■ av, for any a > 0. The above parametrization is therefore suitable for describing the case of the staple links 
extending to infinity. In that case, only the directional information contained in v should enter. For the lattice 
calculations, it is however advantageous to start with an equivalent parametrization in which the structures explicitly 
depend on the staple extent rj and which is still well-defined for v-P = 0. Such a parametrization can be obtained 
from the parametrization above by replacing v ^ -qv and by leaving out the factors (w-P)~". For example, 

^^uiubtr. S, vv, /i) = mjvfli - imle^'-P^P^Kr^VpS.h , (21) 

and analogously for the other Dirac structures. Here we have used lower case amplitudes to distinguish the two 
parametrizations. The relation to the upper case amplitudes is given by 

,2 L D "'^ „„, _ -x {u2 u ID „„, i„„.\-2 



^^{b^b-P,—,j—^,7jvP] ^a,ib',b-P,Tjvb,irjvy,ijvP) 



Bi{b',b-P,^,y^,VV-P] ={vv-Prhib^,b-P,7jvb,{r]vf,vvP), (22) 



vP' [vPf 
v-b 

Vp' {vpy 

where n is the power with which v-P appears in the denominator in front of the corresponding amplitude Bi in the 
parametrization. Notice that the di and bi are functions of all the Lorentz-invariant products of b, P and rjv. For the 
upper case amplitudes, however, we choose to represent the dependence on these invariants in the third and fourth 
argument by 77-independent expressions, in order to facilitate taking the limit rj — ^ ±00. The dependence on the 
Collins-Soper parameter is given by the fourth argument, /(y-Py = — l/(mjvC)^) while the fifth argument, rjv-P, 
characterizes the length of the gauge link and distinguishes between future and past pointing Wilson lines. For the 
calculation of TMDs we work in a frame with 6+ = and = Pt = 0. This leads to a relation that can be expressed 
in Lorentz-invariant form as 



vP 



where 



RiC)^l-^l + C-^ = -^—. (24) 



The relation Eq. ( 23 ) shows that the third argument of the Ai and Bi is not independent of the others in the context 



of TMDs. Moreover, in our lattice calculations, we have to choose the link directions b and v such that Eq. (23) is 



fulfilled. As a side remark, the parameter corresponding to Eq. (23 1 in momentum space is v-k/v-P « x, i.e., the 
amplitudes Ai and Bi acquire an explicit a;-dependence, which has already been pointed out in Refs. pT] |4T|. 
For the F-structures at leading twist, the correlator can be written in the form 

^ ^nnLtr.^^2B+imNe,,hS,AuB (25) 



2P+ 

1 



2P^ 

1 

2P- 



■*u„JbL. = -AA,B+t{{b-P)A-mNibT-ST)} AjB (26) 
-^unsubtr. = imNe.jbj AiB - Aqb - iruNkbi Awe + wjv {(&--P)A - TOjv(bT-S'T)} b, Aub (27) 



where the indices i,j correspond to transverse directions, i,j G {1, 2} (cf. appendix |A] for further details on notation), 
and where we have introduced the following abbreviations for combinations of amplitudes: 

A2B = A2 + RiC^)Bi 

A^B = l4-i?(C')^3 

Mb = l6+(l-i?(C')){Bii+i?(C')Si4} 
AjB = Ar + R{C^)Bi3 
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Aqb = Ag+R{C^)Bi5 

AloB = Alo - (l - i?(C')) {Bi8 - R{C^)B2o} 

AuB = 111 - i?(C')Si7 

Ai2B = Il2-i?(C')^8 (28) 

For later convenience we also define 

AgBm = Aqb " ^mjfb^AnB ■ (29) 

Performing the Fourier transformation and comparing with the decomposition Eqs. Q-Q, we can express the TMDs 
in terms of Fourier-transforms of the above amplitudes. Using the combined amplitudes Aig, the results are of the 
same form as in the straight-link case of Ref . [T3] , 



/i(a;,fc^;C,...,r?i;-P) = 2^ , 
5i(a;,4;C,...,77i>-P) = -2^ Igs + 29, ^ 1 



f ~ 

giT{x,k'^;C,...,'i]v-P) =4m%dk2^ / Atb , 

J J- 

hi{x,k'^-X, . . . ,r]v-P) = -2 J Aqb-hi , 
hidx, C, . . .,-qvP) = 4TO^9fe2, (^j AioB + J 
hir{x,k^T;C---.VvP)^»rn%(^dk2^y Ins, (30) 
except that the abbreviation Jj- is now applied to the w-dependent amplitudes and includes the soft factor: 

I.^ / l^e-"-..^^ / ^,^^(''■P)I^(_6^,fe.p,(5.p)i^(C2)/^2^,_l/(^,^C)^^-P) 
J (271-)^ S{b^;...)J (27r) 



J" 



d{~9) MV^\kT\ 



2(2^) 5(52;...) y (27r) 
Also, there are two further TMDs that are not present in the straight-link case, the T-odd distributions 

/i^(a;,fcT;C, •■•,'7W-P) =4m^afc2 / 1i2b , 

J J- 

/i^(a;,4;C,...,r;«-P) = -4m2^afc2 /" Mb- (32) 

Again the dots ". . ." indicate further parameters that specify the geometry of the soft factor. The T-even distributions 
/i, gi, /ii, giT, and h\j, fuMl 

(x, 4; C, . . . , nv-P) = /^-''^"(x, fc^; C, . . . , -r/^^•P) (33) 
while the T-odd distributions, i.e., at leading twist the Sivers function and the Boer-Mulders function h\; ^ fulfill 

f-°^\x, ; C, . . . , ^V■P) = ~f-°^\x, ; C, . . . , -^VP) (34) 

As a result, T-odd distributions must vanish for 77 = 0, which corresponds to straight gauge links. TMDs for SIDIS 
and DY are obtained for r\v-P — > 00 and r\v-P — > —00, respectively. In the following, we choose v-P > 0, such that the 
SIDIS and DY limits for space-like v can also be written as ri\v\ — > 00 and ?7|w| — > —00, respectively. Equations (30) 
and ([32]) show that certain x-integrated TMDs in Fourier space directly correspond to the amplitudes Ais evaluated 
at b-P = : 

/T'^°^ (b^; C, ■ • ■ , VvP) = 2 l2B(-bT, 0, 0, -l/(mjvC)', VvP)/S{b'; . . .) , 
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~[1](0) 

9i 






.,r]vP) = 


-2 l6s(-b^, 0, 0, -l/(mArC)2, r;t-P)/5(52; . . .) 


9lT 




c,. 


. , rjv-P) ~ 


-2 l7s(-b^, 0, 0, -l/{mNC)^VvP)/§ib''; . . .) 






c,. 


. , rjv-P) = 


-2 l9B,„(-62 , 0, 0, -l/(m^^C)^ . . 


iMm) 


(bx; 


C,- 


. , rjv-P) = 


-2 IioB(-bT, 0, 0, -l/(m^^C)^ 77'^-F)/5(62; . . . 


"it 




C,- 


. , rjv-P) = 


4lni3(-bT, 0, 0, ~l/imNC)^VV-P)/Sib^; ■ ■ ■) , 


J IT 






. , rjv-P) = 


-2 Ii2B(-bT> 0, 0, -l/(mwC)', . . . 




(bx; 


c,. 


. , JJV-P) ~ 


2 l4B(-b|, 0, 0, -l/(TOArC)', rywP)/5(62; . . .) . 



The (derivatives of) Fourier-transformed TMDs f'f\ g[^\ g[l^, \ hfl^^\ hf-^^\ Ivr^^ ^"^^ hf^^'' are naturally 
accessible from the Fourier-transformed cross section of, e.g., SIDIS [33 , and naturally appear in evolution equations, 
see, e.g. 



G. Generalized shifts from amplitudes 



In section |IID| we have given an example that ratios of certain fcT-nioments of TMDs have interesting physical 
interpretations. These ratios, and their counterparts generalized to non-zero bx, are also advantageous from a the- 



oretical point of view: Obviously, the soft factor S cancels in any ratio formed from the objects in Eq. (351, along 
with any F- independent multiplicative renormalization factor jllH131 [?3 . 

Eq. ( 35 1 identifies 2:-integrated derivatives of Fourier-transformed TMDs with simple linear combinations of am- 



plitudes Ai and Bi evaluated at the same values of brp, h ■ P, C, and rjv-P. Forming ratios of these objects thus 

~ fri 

just amounts to taking ratios of linear combinations of the fundamental correlators ^unguj^tr evaluated at the same 
point, i.e., with the same values for 6, P and rjv. For a discussion of the renormalization properties of ratios of the 
objects in (35 1 it is thus sufficient to understand the renormalization properties of (ratios formed from) the correlators 

^Ifiubtr. = 1 (^.^1 m^Uq{h) \P,S). 

Analytical studies of the operator g(0) TUq{b) in the continuum [421 - 147] suggest that for 6^ 7^ the renormalization 
factors are multiplicative and F-independent. The basic reason is that the quark field operators are at different 
locations and undergo wave function renormalization separately. We will assume here that our lattice representation 
of q{0)rUq{b) is renormalized multiplicatively independent of F as long as we keep larger than a few lattice 
spacings. A more detailed discussion and numerical studies of the renormalization properties of this operator can 
be found in Ref. |13j . It remains an interesting task for the future to perform a more thorough treatment of 
non-local operators on the lattice. Under the assumption of multiplicative renormalization, generalized shifts such 

as {ky)Tu 

= mAr/j!jl^'^^''//j^'^°^ can only depend on b'^, ( and on the staple extent rjv-P. AU other 
renormalization and soft factor related dependences cancel out in the ratio. In this work, we will present numerical 
results for the following generalized shifts: 



■P) 



...) 



{h^)Tij\b^-C,r\v-P) = mjv t,^.,^. — -111^— ^ 

/W(°)(b^;C,...,^«-P) AM-bl,Q,0,-l/{m^Cr,VV-P) 

bl=o J dxj (fk^ky <P^''^\x,kT,P,S; 
^ JdxJd^kT <i>^'r+]{x,kT,P,S; 

ht^'^^'\b^^;C...,VvP) AiB{-b^T,0,0:-^/{rnNC)^Vv 

(6^; C, . . . , Vv-P) """^ A2B{-b^T, 0, 0, -1/ {mNCY ,VvP: 



St = (1,0) 



(36) 



iky)uAb,;C,VvP) = --^p)(,2^^^ 



Ea=±1 JdxJ d^kT $[7 + +.--^ + 7=] {x, kT, P,S-.. .) 



St 



= (1,0) 



(37) 
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(..)..(.,;C,...P) ^ --;[.Ko)(,.^^...^^„.^)--^^^(_,.,o,0,-l/(™.0^..^ 
b|=o / / d^fcT fex ^I-^^+^-^^-^'l (a;, kT,P,S;.. .) 

^ JdxJd^kT $[7+ + A7+7'^](2;,feT,P,S';...) 



/Tl^°)(b2 ; . . . , ,7t;.P) A2B(-bT, 0, 0, ~l/(mwC)^ ^^'■^) 



b; 



^ JdxJd'^kT $['>'+l(a;,fcT,P,5;...; 



St = (1,0), A = l 



St = (1,0), ST = (1,0) 



(38) 



(39) 



• The "generalized Sivers shift" (fcj^)^'™'''' — {ky)Tu has aheady been discussed in section II D It is T-odd, i.e., we 
expect to obtain results of opposite sign in the SIDIS and DY limits r]v-P — ^ oo and r]v-P — >■ —oo, respectively. 
The generalized Sivers shift describes a feature of the transverse momentum distribution of (unpolarized) quarks 
in a transversely polarized proton. In the formal limit ^t = it measures the dipole moment of that distribution 
orthogonal to the polarization of the proton. 

• The "generalized Boer-Mulders shift" (fcj,)^'^ — {ky)uT is also T-odd and addresses the distribution of trans- 
versely polarized quarks in an unpolarized proton. In the limit = 0, the Boer-Mulders shift describes the 
dipole moment of that distribution orthogonal to the polarization of the quarks. Note that we use a sum over 
proton helicities X]a=±i ( |37P to represent the unpolarized target nucleon. 

• The generalized shift (kj.)^^'^ — {kx)TL attributed to the "worm gear" function gix quantifies a dipole deforma- 
tion of the transverse momentum distribution induced by the correlation of the quark helicity and the transverse 
proton spin. Unlike the Sivers and the Boer-Mulders shifts, it is a T-even quantity, i.e., the SIDIS and DY limits 
■qvP — T' ±oo are expected to be the same. This shift has already been studied in lattice QCD using straight 
gauge links [TTHI3]- We are interested to see by how much this "process independent" result obtained at = 
differs from the results calculated with SIDIS- and DY-type gauge links in the limit rjv-P zLoo. 

• The ratio /if '^°V/f ''"^ 

can be identified with a "generalized tensor charge" . Clearly, it is also a T-even quantity, 
i.e., no differences are expected between the SIDIS and DY limits rjv-P — > ±oo. We have studied already in 

[TTHT5] on the lattice using straight gauge links. As doesn't involve any fc-weighting and is directly 

related to the well-known transversity and unpolarized distribution functions, we expect it to be a particularly 
clean observable. It therefore qualifies as a very good candidate for our study of the 77|w|-dependence of T-even 
observables, in particular the transition from straight to staple-shaped gauge links. 

The framework laid out above provides the basis for our numerical lattice calculations described in the next section. 
Before proceeding, it is worth reiterating the logic underlying our approach. Recognizing that the generic kinematics 
for which TMDs are defined are space-like, with light-like separations representing a special limiting case, we proceed 
by considering kinematics off the light cone from the start. We again emphasize that, whereas we thoroughly examine 
the behavior of our data as the kinematics are pushed in the direction of the light cone, statements about formal 
properties of the light-cone limit lie beyond the purview of this investigation. Having parametrized the relevant 
nonlocal matrix element in terms of Lorentz-invariant amplitudes, cf. section ["nF[ we choose to perform its evaluation 
in a Lorentz frame in which the operator under consideration is defined at one fixed time. There is no obstacle to 
this choice in view of the space-like separations entering the original definition of the matrix element. In this frame, 
we cast the computation of the matrix element in terms of a Euclidean path integral, which we evaluate employing 
lattice QCD, as detailed in the next section. 



III. LATTICE CALCULATIONS 
A. Simulation setup and parameters 

The methodology we use to calculate the non-local correlators on the lattice has been described in detail in Ref . [13] , 
except that we now extend this method to staple-shaped links. Again, we employ MILC lattices [48l|49] that have been 
previously used by the LHP collaboration for GPD calculations |50j : however, compared to our previous work with 
straight gauge links, we now go to lighter pion masses and make use of the coherent proton and anti-proton sequential 
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rhu^d rhs 


xf 




a (fm) 


m?*P (MeV) 


niff (GeV) 


T^conf. 


#meas. 


0.01 0.05 
0.01 0.05 
0.02 0.05 


28^ X 64 
20^ X 64 
20^ X 64 


6.76 
6.76 
6.79 


0.11967(14)(99) 
0.11967(14)(99) 
0.11849(14)(99) 


369.0(09)(35) 
369.0(09)(35) 
518.4(07)(49) 


1.197(09)(12) 
1.197(09)(12) 
1.348(09)(13) 


273 
658 
486 


2184 
5264 
3888 



TABLE L Lattice parameters of the n/ — 2+1 MILC gauge configurations [481 149j used in this work. The lattice spacing a 
has been obtained from the "smoothed" values for ri/a given in Ref. [52] and the value ri = 0.3133(26) fm from the analysis 
of Ref. [53]. The first error estimates statistical errors in ri/a, the second error originates from the uncertainty about ri in 
physical units. We also list the pion and the nucleon masses determined in Ref. [5T] with the LHPC propagators using domain 
wall valence fermions. The first error is statistical, the second error comes from the conversion to physical units using a as 
quoted in the table. Note that the masses quoted here in physical units differ slightly from those listed in Refs. |50l I51j . 
because these references use a different scheme to fix the lattice spacing. The second to last column lists the number of gauge 
configurations and the last column shows the resulting number of measurements for the calculation of three-point functions 
achieved by means of multiple locations for source and sink. 



propagators of Ref. [5T] to increase our statistics. The new LHPC data set offers forward propagators at four different 
source locations on each gauge configuration. Moreover, coherent proton and antiproton sequential propagators have 
been calculated, each one implementing simultaneously four nucleon sink locations per gauge configuration. This way 
it is possible to conduct eight measurements of a three-point function on each gauge configuration in well separated 
areas of the lattice, boosting statistics significantly. The source-sink separation has been chosen to be nine lattice 
units. The simulation parameters are summarized in Table [l] 



B. Nucleon momenta, choice of link directions, and extraction of amplitudes 

For all of the ensembles listed in Table |l| nucleon momenta P — and P — 27r/(ai) • (—1,0,0), implemented via 
corresponding momentum projections in the sequential propagators, were available. In addition, sequential propaga- 
tors were produced corresponding to the nucleon momenta P — 27r/(aL) • (—2, 0, 0) and P = 2tt/ {aL) •(1,-1,0) for the 

mu,d = 0.02 ensemble only. We extracted the matrix element 5|fJ,ubtr. '5';Cfc) = ^ (P, S\ q{0) T U[Cb] g(6) |P, S) 
from plateaux in standard three-point function to two-point function ratios, for a complete basis of F structures and 
nucleon states polarized in the 3-direction. The nucleon momenta P, quark separations b and corresponding staple- 
shaped gauge link paths Cf, used on the lattice in the present investigation are listed in Table |lTj The link path Cf, is 
characterized by the quark separation vector b and the staple vector rjv, cf. Fig. [T] The range of 77 studied was always 
chosen to extend from zero to well beyond the point where a numerical signal ceases to be discernible. Furthermore, 
it should be noted that in the case of either b or v extending into a direction in a lattice plane which forms an angle of 
7r/4 with the lattice axes spanning the plane, there are two optimal approximations of the corresponding continuum 
path by a lattice link path; e.g., if one denotes the lattice link vector in i-direction as e^, then b = 2(ei + 62) is 
equally well approximated by the sequence of links (61,62,61,62) as by the sequence (62,61,62,61). As far as b is 
concerned, in such a situation, our calculations always included both optimal link paths. However, in the case of v, 
in these situations, only one of the two link paths was included. To be specific, in the instances of r]v = ±n'{ei ± 6^) 
quoted in Table |Tlj the link path always departs from the quark locations in i-direction, not 1-direction. This is a 
shortcoming of the discretization which breaks the manifest T-transformation properties present for the continuum 
staple; presumably it is responsible for the problematic mixing of T-even and T-odd amplitudes which we observe in 
our analysis in the case of staple directions off the lattice axes. While we expect a symmetry-improved calculation 
including both optimal link paths to avoid this issue, with the presently available data, we find that we need to impose 
explicitly T-odd/T-even symmetry in the system of equations from which we extract the amplitudes whenever v does 
not coincide with a lattice axis. 

— fr] 

In practice, the overdetermined system of equations which we solve in order to relate the matrix elements 'I'unsuMr 



to t he corresponding amplitudes is set up in terms of the quantities ai,bi, cf. Eq. (22 1 in conjunction with Eqs. (16) 



(20 1. This form is suited to include the case C = Oj where the sign of the prefactor in front of bi depends on whether 



the limit rjv-P = is approached from the SIDIS or the DY side. 
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b/a 


T]v/a 


■ aL/[2-K) 


Notes 


n ■ (0,0, l),n = -7, ... ,7 


±n' ■ (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


rhu,d = 0.02 ensemble only 


±n' • (1, 1,0) 


(-1,0,0) 




(-2,0,0) 


'rhu,d = 0.02 ensemble only 


±n' • (1,0,0) 


(1,-1,0) 


'm,u,d = 0.02 ensemble only 


±n'- (1,-1,0) 


(1,-1,0) 


= 0.02 ensemble only 


n- (0, l,Q),n = -7,..., 7 


±n' • (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


Thu,d = 0.02 ensemble only 


±n' ■ (0,0,1) 


(-1,0,0) 




(-2,0,0) 


rhu,d = 0.02 ensemble only 


±n' • (1,0, 1) 


(-1,0,0) 




(-2,0,0) 


'rhu,d = 0.02 ensemble only 


n ■ (0, 1,1), n = -2, ... ,2 


±n' • (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


WT'u.d = 0.02 ensemble only 


n- (0,-l,l),n = -2,...,2 


±n' • (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


Thu,d = 0.02 ensemble only 


±(0,3, ±2) 


±n' • (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


'mu,d = 0.02 ensemble only 


±(0,4, ±2) 


±n' • (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


= 0.02 ensemble only 


±(0,4, ±3) 


±n' • (1,0,0) 


(0,0,0) 




(-1,0,0) 




(-2,0,0) 


'rhu,d = 0.02 ensemble only 


n ■ (1, l,0),n = -4, . . . ,4 


±n'- (1,-1,0) 


(1,-1,0) 


Thu,d = 0.02 ensemble only 



TABLE II: Sets of staple-shaped gauge link paths and nucleon momenta P used on the lattice. Gauge link paths are character- 
ized by the quark separation vector b and the staple vector rjv, cf. Fig.[l] The surveyed range of rj, parameterized in the table by 
the integer n', was always chosen to extend from zero to well beyond the point where a numerical signal ceases to be discernible. 
The maximal magnitude of the CoUins-Soper parameter C, attained in these sets is |^| = 0.78, for P-aL/{2-K) = (—2, 0, 0) paired 
with riv/a = ±n' ■ (1,0,0). 

C. Numerical Results 

1. The generalized Sivers shift 

In the following, we concentrate on results for the isovector, u — d quark combination, because in this case contribu- 
tions from disconnected diagrams and possible vacuum expectation values cancel out. The errors shown are statistical 
only. At the present level of accuracy in this exploratory study, we set aside a quantitative analysis of systematic 
errors. We use the central values for the lattice spacing a as given in Table |T] to convert to physical units. For m^r, we 
consistently substitute the value of the nucleon mass as determined on the lattice, rather than the physical nucleon 
mass. 

Figures [s] to It] show our results for the generalized Sivers shift, {ky)^^'^^. We begin with a discussion of its 
dependence on fne staple orientation, i.e., SIDIS- or DY-likc, and the staple extent, ri\v\, as displayed in Fig. [s] for a 
Collins-Soper evolution parameter of (■ = 0.39 and a pion mass of = 518 MeV. As mentioned before, the T-odd 
Sivers function must vanish for 77|?;| = 0, i.e., a straight Wilson line between the quark fields, but non-vanishing 
results are allowed (and generally expected) for non-zero staple extents. Furthermore, the T-odd observables are 
anti-symmetric in ri\v\, so we expect the Sivers shift to be of the same size but opposite in sign for the SIDIS and 
the DY cases. This is exactly what we find in, e.g., Fig. |3a[ showing the shift for a quark-antiquark distance of a 
single lattice spacing, |6t| = la- The aforementioned features are realized in form of a curve that is reminiscent 
of a hyperbolic tangent. We stress that the observed zero crossing with a change in sign is directly caused by the 
underlying gauge-invariant operators and their symmetry properties, and hence represents a consistency check of our 
calculation rather than any sort of a prediction. 

Remarkably, already as \ri\\v\ approaches values of ~ 6a, we find that the Sivers shift stabilizes and reaches specific 
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FIG. 3: Extraction of the generalized Sivers shift on the lattice with m,r = 518 MeV using a lattice nucleon momentum 
|piat| _ 2-K/{aL) ~ SOOMeV at the corresponding maximal Collins-Soper evolution parameter (J = 0.39. The continuous 
horizontal lines are obtained from two independent averages of the data points with staple extents in the ranges ri\v\ — 7a. .12a 
and ri\v\ = —12a.. — 7a, respectively. The outer data points shown with empty symbols have been obtained from an anti- 
symmetrized mean value of these averages, i.e., the expected T-odd behavior of the Sivers shift has been put in explicitly. 
These outer data points are our estimates for the asymptotic values at ?7|ii| — >■ ±oo and thus represent the generalized Sivers 
shifts for SIDIS and DY. Error bars show statistical uncertainties only. Figures |(a)| and |(b)| have been obtained with rather 
small quark field separations |bT| = Iffl and 2a. Therefore, they might be affected by significant lattice cutoff effects. 



plateau values. Apart from finite volume effects, in particular wrap-around effects due to the periodic boundary 
conditions on the lattice, we see no reason to expect that once a plateau has been reached, the value of the shift 
would significantly change as |77||w| — ^ c». To obtain first estimates for staple-shaped Wilson lines that have an 
infinite extent in u-direction, we therefore choose to average the shifts in the plateau regions \ri\\v\ = 7a . . . 12a, as 
illustrated by the straight lines. Clearly, as Ibx] increases from 0.12 fm in Fig. 3a to 0.47 fm in Fig. 3d the signal- 
to- noise ratio decreases as we approach larger values of |?7||w|. For smaller |ry| \v\ the statistical uncertainties are 
much smaller, and the corresponding values tend to dominate the averages when the errors are taken into account as 
weights. At the same time, however, these statistically dominating data points are more likely to introduce systematic 
uncertainties related to the (unknown) onset of the "true" plateau region and the corresponding starting value for 
the averaging procedure. Therefore, in order to avoid a too strong bias from the data at smaller |f7||w|, we do not use 
the respective statistical errors as weights in the averaging. Our final estimates for the Sivers shift are obtained from 
the mean value of the SIDIS and DY averages and by imposing antisymmetry in r]\v\. The results are displayed as 
open diamonds at r]\v\ = ±oo in Fig. [3j The dependence of these results on |6t| is shown in Fig. |4] In summary, 
for (■ = 0.39 and I^tI — 0.12 . . . 0.47 fm, we find a sizeable negative Sivers shift for w — d quarks in the range of 
{kyfu-T'^'''''^ = -0-3 ... - 0.15 GeV. 

Next, we turn to the dependence of our results on the Collins-Soper evolution parameter ^. In Fig.[5j we consider two 
"extreme" cases, namely, a vanishing ^ as well as the largest ( 0.78 that we could access in this study. While we find 
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FIG. 4: Generalized Sivers shift as a function of the quark separation Ibxl for the SIDIS case — oo), extracted on the 

lattice with = 518MeV for = 0.39. The data points lying in the shaded area below |bTl ~ 0.25 fm might be affected by 
significant lattice cutoff effects. Error bars show statistical uncertainties only. 
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FIG. 5: Generalized Sivers shift on the lattice with — 518 MeV for a quark separation of three lattice spacings^Jbr j = 3a 



0.36 fm, extracted at = and at our highest value of the CoUins-Soper evolution parameter, = 0.78. Figure (b) has been 
obtained from nucleons with momentum IP''^'] = 2 x 2iT/{aL) ~ 1 GeV on the lattice. Error bars show statistical uncertainties 
only. 



rather precise values for the Sivers shift ioi ( — with a well-defined plateau^ for \t]\\v\ > 6a in Fig. 
and uncertainties quickly increase with |77||w| for <^ = 0.78 in Fig. 



of a plateau on the right hand (SIDIS) side of Fig. 5b Following t 
find that the estimated values at |?7||w 
Figure 



5b 



5a 



fluctuations 



In particular, it is difficult to identify the onset 
le averaging procedure described above, we however 



6a 



c» for the two extreme cases of C agree within uncertainties, 
shows the Sivers shift as a function of (, for |6t| = 0.36 fm and a pion mass of — 518 MeV. Within 
the present uncertainties, we observe a statistically significant negative shift; however, it is not possible to identify a 
clear trend of the data points as C increases. With respect to data points obtained for staple link directions v off the 
lattice axes, i.e., ( w 0.55 in Fig. 



6a 



we note again that we need to impose the T-odd/T-even (anti-)symmetry already 
when we solve our system of equations, in order to avoid problematic mixings of T-even and T-odd amplitudes. As 
already mentioned further above, we expect this to become unnecessary in the case that lattice symmetry improved 
operators (see Appendix D of Ref. [13]) are used. 



We find it very interesting to note that the contribution from A12 alone in the numerator of Eq. (36) (rather than 



^ Note that, in the case at hand, C = corresponds to P = 0, so that one cannot identify a "forward" or "backward" direction. Hence, 
there is only a single branch in r]\v\, the sign of which is a matter of definition, see also the discussion further below in the text. 
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FIG. 6: Generalize d Si vers shift on the lattice with — 518 MeV for a quark separation of three lattice spacings, |bTl = 3a = 
0.36 fm. In Figure (a) we show the ("-dependence of the generalized Sivers shift, depicting both the full result and the result 
obtained wit h ju st A12 i n th e numerator. The data points correspond to those displayed in the SIDIS limit 77|?;| — >■ cxa in plots 
such as Fig. (b) Figure (b) shows the 77-dependence at — 0.39 for both the full result (diamonds) and the contribution from 
amplitude A12 in the numerator (squares). Asymptotic results corresponding to SIDIS and DY have been extracted as in Fig. 
[3] Error bars show statistical uncertainties only. 



Ai2b), illustrated by the open squares, is essentially compatible with zero within errors for all accessible values of C. 
The main contribution to the transverse shift therefore comes from — i?(C^)i3g = —ri{v ■ P) R{(^)bg, (see Eqs. (22)), 
i.e., the amplitude &g. Note again that, on the lattice, we employ expressions in terms of the lower-case di and bi 
amplitudes, e.g. bs, as they are well defined even when ( 0. In this limit, v- P ^ 0, and hence the prefactor behaves 
as —{v ■ P) R{C^) arriM- The sign of the prefactor of bi depends on whether one approaches the limit u • P — !• 
from the SIDIS or the DY side. 



An example that explicitly shows the relative smallness of A12 is given in Fig. 6b for C = 0.39 and |5t| = 0.36 fm. 
While A12 as a function of r/|i;| shows the typical behavior expected for a T-odd amplitude, it represents only about 
10% of the total contribution for, e.g., |?7||w| = 6a. 

As one of our central results, we show in Fig. [rjthe Sivers shift as a function of C, for all considered ensembles, as 
before for a fixed I^tI = 0.36 fm. Within statistical uncertainties, the data points for the two different pion masses 
m^r = 369 MeV and ttItt — 518 MeV, as well as the spatial lattice volumes V « (2.4 fm)"^ and V « (3.4 fm)^, are overall 
well compatible. Apart from the less well determined data point at C 0.55, we find a clearly non-zero negative Sivers 
shift in the range {ky)^^^^'^^^^^ = —0.48 ... — 0.2 GeV. Together with the relatively mild bx-dependence at smaller 
bx, cf. Fig. [4] this provides strong evidence that the (x- and fcx-moment of the) Sivers function considered here 
is sizeable and negative for u ~ d quarks. Our preliminary separate data for u- and for d-quarks (not shown in this 
work) furthermore indicate that f-^rj.^ < and f^r^ > 0. Although our results for the T-odd Sivers effect are still 
subject to many systematic effects and uncertainties, it is interesting to note that they are overall well compatible with 
results from a phenomenological analysis of SIDIS data [541 155| , as well as arguments based on the chromodynamic 
lensing mechanism by Burkardt [36l[56lE7|. It should also be noted that, in a recent twist-3 analysis of single spin 
asymmetries from RHIC experiments, a possible discrepancy has been found with respect to the signs [5H] . 

We stress again that fully quantitative predictions for, or a comparison with, phenomenological and experimental 
TMD studies employing QCD factorization would require lattice data for much larger CoUins-Soper parameters, 
(■ ^ 1. With = — (w-P)^/(f ^m^j, the limit ( —>■ 00 corresponds to the limit of a light-like staple direction v, or an 
infinite rapidity yy ^ —00 in Eq. (15). For the shifts and ratios defined in Eqs. (36)- (39 1, where the soft factors in the 
TMD definitions [HEllTj cancel out, large values of ( can also be accessed through large nucleon momenta. Clearly, 
the limit of an infinite Collins-Soper parameter is in practice not accessible on the lattice, so that we have to rely on 
results for a limited range of (, as for example in Fig. [t] From the perturbative prediction for the C dependence, cf., 
e.g., Ref. [6], we would expect that the ratios of TMDs should become independent of ^ as (" — >■ cxd. It would be very 
interesting to investigate this on the basis of future lattice results for larger hadron momenta and with substantially 
improved statistics. 
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FIG. 7: Comparison of the (["-evolution of the generaHzed Sivers shift at |?>t| = 3a = 0.36 fm for the three different lattices 
listed in Table [T] Filled symbols correspond to the full SIDIS result. The data points with open symbols have been obtained 
with only A12 in the numerator. Error bars show statistical uncertainties only. 



2. The generalized Boer-Mulders shift 



We now turn to the second pr omi nent T-odd TMD, the Bo er-M u lder s function. Our results for the generaUzed 
Boer-Mulders shift (ky)^^^ (Eq. (37)) are summarized in Figs, 
is shown in Fig. 



to 



8d 



A typical example for the ry|u|-dependence 



for a pion mass of = 518 MeV, C = 0.39, and |6t| ~ 0.36 fm. Apart from the magnitude of 
the shift, the results are very similar to what we have found for the Sivers shift in Fig. [3c] above, with indications 
for plateaus for > 6a. Figure 8b illustrates the dependence on |6x| for the SIDIS case. Although the central 



values indicate some trend towards values smaller in magnitude as [bxl increases, the somewhat large uncertainties 
and fluctuations at larger |6t| prevent us from drawing any strong conclusions. In the range of Ibxl ~ 0. . .0.4fm, 
we find a clearly non-zero negative Boer-Mulders shift of {ky)^^^^^^^ ~ —0.17 ... — 0.1 GeV, for ^ = 0.39 and the 
The ^-dependence for m.^- 



8c 



As for the Sivers 



given pion mass, ihe C,-aepenaence tor m.^- = 518MeV and Ibxl = 0.36 fm is shown in Fig. 
shift, it is interesting to note that the contribution from the A amplitude, in this case A4, given by the open squares, 
is mostly compatible with zero within errors, while the main signal is coming from —R{Q^)B^ — —7] (v ■ P) R(C,^)b^, 
cf. Eqs.'([28l). 



Finally, a comparison of the results and their ^-dependences for the three different lattice ensembles is provided 



in Fig. 8d for a fixed |bx| — 0.36 fm. We find that most of the data points for the two pion masses and the two 
volumes are well compatible within uncertainties, with central values of {ky)^^^^^^ ~ —0.2... — 0.1 GeV. While 

the central values show little dependence on C, the errors have to be significantly reduced before any extrapolations 
towards a large CoUins-Soper parameter may be attempted. In summary, for the given ranges of parameters, our 
results indicate that the Boer-Mulders function is sizeable and negative for u~ d quarks. Our data for the individual 
u- and d-quark contributions (not shown) furthermore indicate that /i^'" < and h^'^^ < 0. Interestingly, these 
preliminary results are well compatible with a recent phenomenological study of the Boer-Mulders effect in SIDIS 
|59j . as well as an earlier lattice QCD study of tensor generalized parton distributions |60| in combination with the 
chromodynamic lensing mechanism [61j . 



3. T-even TMDs: The transversity hi 



In the previous sections, we have discussed the T-odd Sivers and Boer-Mulders distributions, in particular their 
emergence in the transition from straight to staple-shaped gauge links, i.e., as r]\v\ changes from zero to large positive 
or negative values. A natural question to ask is, what is the infiuence of final state interactions, which we mimick 
on the lattice with the staple-shaped links, and which are essential for the appearance of T-odd distributions, on the 
T-even TMDs? More specifically, we would like to see whether and how the T-even distributions, which are generically 
non-vanishing already for straight gauge links, change during the transition to finite staple extents. This is also of 
considerable interest with respect to the much less involved lattice studies of (T-even) TMDs using straight gauge 
links that we have presented in jl3j . As we will show, there is only little difference in the transition to staple-shaped 
links, such that our previous results might be of greater phenomenological importance than initially expected for the 
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FIG. 8: Generalized Boer-Mulders shift. 

(a) r;j«|-dependence at — 518 MeV for ( — 0.39, I^t] ~ 3a — 0.36 fm. 
Asymptotic results corresponding to SIDIS and DY have been extracted as in Fig. [s] 

(b) ]|bT|-dependence of the SIDIS results at m„ = 518 MeV, ( = 0.39. 

(c) fc-dependence of the SIDIS resuhs at m„ = 518 MeV, I^tI = 0.36 fm. 
Empty squares correspond to the ratio with A4 in the numerator only. 

(d) Comparison of the (^-dependence of the SIDIS results obtained from the three different lattice ensembles listed in Table |l] 
All error bars show statistical uncertainties only. 



straight "process-independent" gauge link structures. 

A suitable observable for investigating these questions is the "generalized tensor charge" given by the ratio of the 

defined in Eq. The ri\v\- 

for different Ibx] of 0.12, 0.24, and 0.36 fm, a pion 



(lowest a;- moments of the) transversity to the unpolarized distribution, h}^^^^^ / f\^^^°\ 



dependence of this transversity ratio is displayed in Figs. 9a to 9c 



the full range of accessible |?7||t;| in Figs. 9a and 9b 



mass of m^r — 518 MeV, and ^ = 0.39. We find it quite remarkable to see that stays nearly constant over 

within comparatively small statistical errors. For Ibxl = 0.36 fm, 
?7||t;| where the signal-to-noise ratio quickly decreases. In all cases, 
8a. As in the previous sections, we choose to average over the data 
±00 from the mean of the DY and SIDIS averages 



we see little dependence apart from larger values of 
we find indications for plateaus from |77||w| ~ 3a . 
in the plateau regions (solid hues), and obtain estimates for \r]\\v\ 
by imposing the symmetry condition in r]\v\. The corresponding results are illustrated by the open diamonds. In all 
considered cases, differences between |77||u| — > ±00 and = are barely visible within uncertainties. In other 

words, lattice data for simple straight gauge links provide already a very good estimate for the phenomenologically 
interesting case of infinite staple extents, at least in the covered ranges of C and not too large |bT|. 

The Ibxl-dependence of our estimates for the transversity ratio at \ri\\v\ = ±00 is displayed in Fig 9d for m^r — 
518 MeV and C = 0.39. We observe a small, approximately linear rise of in total about 20% as |6t| increases 
from 0.12 fm to about 0.6 fm. This is in agreement with our previous observation of a flatter |bT|-dependence of 
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FIG. 9: 



The dependence of the transversity ratio Eq. (391, on the staple extent ri\v\, obtained at rrin = 

518 MeV, ( — 0.39 for three different quark separations |bTj — la — 0.12 fm, 2a — 0.24 fm and 3a — 0.36 fm. Asymptotic results 
corresponding to SIDIS and DY have been extracted as in Fig. [S] except that we assume an even behavior of hi to obtain 
the data points plotted as empty symbols at ri\v\ — > ±oo. The averages (lines) are obtained from the data points with staple 
extents in the ranges i]\v\ = 3a. .8a and 77|i;| = —8a.. — 3a, respectively. Figure (a) might be affected by significant lattice cutoff 
effects due to the small quark separation |6t| = a. 



(d) The transversity ratio ^j^""'//]'^'^"^ as a function of the quark separation |bT| from the SIDIS results extracted on the 
lattice with = 518MeV for — 0.39. The data points lying in the shaded area below |6t1 ~ 0.25 fm might be affected by 
lattice cutofi' effects. Error bars show statistical uncertainties only. 



the amplitude Ag^"^ compared to in Ref. [T3] on the basis of straight gauge links'^. Remarkably, a naive 

linear extrapolation of the data to |6t| — would give a value for the tensor charge, g^^'^ = Jdx (Pk^ hi{x,k^) — 
/iW(°'(bT = 0), of ^^""^ « 1.1, which is in very good agreement with the direct lattice calculation of this quantity using 
a renormalized local operator that has been presented in for the same ensemble, for a scale of /i^ = 4 GeV^ in the 
MS scheme''. 

Figure 10a shows the transversity ratio as a function of the CoUins-Soper parameter, for the same pion mass as 
before but a fixed 1^x1= 0.36 fm. The C-dependence turns out to be rather flat over the full range of accessible values. 
It is interesting to note that, in contrast to the T-odd distributions discussed before, th e amplitude Ag„i (open circles) 
provides ~ 100% of the total results, while the contribution from i?(C^)Bi5, Eq. (28 1, as well as from Bn through 
Eq. (29), is negligible within errors, over the full range of C- 



^ Note again that b in the present work corresponds to — / in |13j . 
* Note that the tensor charge is denoted by (l)iq in I62| . 



21 



1.5 



I 1.0 

1 



0.5 



0.0 



0.0 



(a) 



J total 

^ contrib. Ag 



hi, u-d - quarks 

0.2 0.4 0.6 0.8 1.0 



1.2 



Ibj-I = 0.36 fm 
m„ = 518 MeV 



1.4 



1.5 



^ 1.0 



0.5 



0.0 



0.0 



FIG. 



10: Evolution with respect to ^ for the transversity ratio /i j^"-"' //j^"-"^ 



0.2 



0.4 



(b) 



^ m;,=518MeV 20^ 
^ m;r=369MeV 20^ 
^ m„=369MeW 28^ 



IbrI = 0.36 fm 



hi, u-d - quarks 

0.6 0.8 1.0 



1.2 



1.4 



at a quark separation of |6t| — 3a = 0.36 fm. 
Figure [TOa| shows the SIDIS results obtained at m,r = 518 MeV. The sohd data points correspond to the full result, and empty 
symbols to the result obtained with just A^rn in the numerator. Figure 
Table n 



10b 



displays the full results for all ensembles listed in 



Finally, we show a comparison of our results for obtained for the different ensembles in Fig. 



10b 



As before, the data points for the two values of the pion mass and the different volumes agree within uncertainties. 
On the basis of the comparatively good signal-to-noisc ratio for this observable, we conclude that the ^-dependence 
is in this case rather flat and very well compatible with a constant behavior, (/if ''°V/P^°^)""''(C) ~ 1.2 ± 0.1, at 
least for < 0.8 and the given parameters. It would be interesting to investigate in future lattice studies whether this 
constant behavior persists as one approaches larger Collins-Soper parameters. 



4- T-even TMDs: The generalized worm gear shift from gir 



As a final example, we stu dy in this section the generalized shift defined in Eq. (38), which is essentially given by 

the T-even TMD gir- Figure 11 shows (sfj^'^V/l^'*"^)"^'^ as a function of f]\v\ for = 518 MeV, C = 0.39 and two 
values of I^tI- Within uncertainties, we observe, as expected, an approximate symmetry with respect to the sign of 
rylvj. Furthermore, we find overall only little dependence on the staple extent for |6t| — 0.12 fm. At larger Ibxl, the 
signal-to-noise ratio quickly decreases as |?7||w| becomes larger. Still, we find indications that the results stabilize in 
the region \r]\\v\ — 3a . . . 8a, which we choose as our plateau region for the computation of average values, as discussed 
in the previous section. 

As before, the averages serve as approximations for the asymptotic results at r]\v\ — ±oo, i.e., corresponding 

Although 
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to infinite staple extents. The dependence of these asymptotic values on |6x| is displayed in Fig. 
a small curvature in the central values can be observed, the results are overall rather stable within errors, with 
(fex)rd = (Mt'V/?'^"^)""' « 0.16 ... 0.21 GeV. 

In Fig. |13a[ we show the dependence of the generalized gix shift on the Collins-Soper parameter, for a pion mass of 
m^r = 518 MeV and a fixed |bx| = 0.36 fm. As C increases, one observes a slight trend towards values that are smaller 
in magnitude, although it is difficult to draw any strong conclusions in view of the present uncertainties. Similar to 
the case of the transversity ratio of the previous section, we find that essentially the full signal is due to the amplitude 
At, while the c ontri bution from i?(^^)i?i3, cf. Eq. (28), is compatible with zero within errors. 

Finally, Fig. 13b gives an overview of our results as functions of ^, obtained for the three considered ensembles, 
for |6t| = 0.36 fm. Apart from C = 0, the data points for the two pion masses and volumes clearly overlap within 
uncertainties. Taking into consideration the results for m^r — 369 MeV and a spatial volume of 28'^ (given by the 
filled diamonds), the data are overall compatible with a constant behavior, although more statistics is necessary to 
establish a clear trend in ^. Altogether, we observe a sizeable positive generalized transverse shift in the range of 

{kx)uld = (Mt^V/?'^"^"'* ~ 0.15... 0.25 GeV, for C = 0...0.8 and the given parameters. We note that these 
values are in good agreement with our previous analyses on the basis of straight gauge links [TH [13] . 
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FIG. 11: Dependence of the generalized gir shift on the staple extent ri\v\, obtained at = 518 MeV, — 0.39 for two 
different quark separations |bT| ~ Id = 0.12 fm and Ibxl = 3a = 0.36 fm. Asymptotic results corresponding to SIDIS and 



DY have been extracted as in Figs. 9a tol9cl Error bars show statistical uncertainties only. Figure (a) might be affected by 



significant lattice cutoff effects due to the small quark separation Ibxj — a. 
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FIG. 12: Generalized giT shift for I?;] \v\ = cxd as a function of the quark separation |bT| from the SIDIS and DY results extracted 
on the lattice with = 518 MeV for = 0.39. The data points lying in the shaded area below |bTl ~ 0.25 fm might be affected 
by lattice cutoff effects. Error bars show statistical uncertainties only. 



IV. SUMMARY AND CONCLUSIONS 



We have presented an exploratory study of quark transverse momentum distributions in the nucleon in full lattice 
QCD employing non-local operators with staple-shaped gauge links (Wilson lines). Compared to our earlier works 
[m [13] , the use of staple-shaped instead of straight link paths allowed us for the first time to systematically access the 
naively time-reversal odd (T-odd) observables, in particular the amplitudes related to the Sivers and the Boer-Mulders 
TMDs. In the framework of QCD factorization theorems, the path dependence corresponds to a process dependence 
that leads to the famous sign difference between the T-odd TMDs for the SIDIS and the DY processes. In our study, 
we were able to distinguish the SIDIS and DY cases through the relative orientation of the nucleon momentum P 
and the vector rjv that characterizes the direction and extent of the staple on the lattice, cf. Fig. |lb[ It is important 
to keep in mind that TMDs defined with non-light-like staple vectors v, as required on the lattice, will additionally 
depend on the Collins-Soper evolution parameter, here denoted by (. In order to avoid additional soft factors in the 
formal definition of the TMDs, we have concentrated on the Sivers, Boer-Mulders, and worm-gear (gir) generalized 
transverse momentum shifts and the generalized tensor charge. Since the generalized shifts and tensor charge are 
defined in terms of ratios of TMDs, potential soft factors as well as the renormalization constants cancel out. Our 
numerical results, obtained for three different ensembles with pion masses TOtt = 369 MeV and m-^ = 518 MeV, as 
well as spatial lattice volumes of 20'^ and 28'^, are very promising: We find clearly non-zero, sizeable signals for all 
observables we considered. The expected anti-symmetry (change of sign) for T-odd quantities in r]\v\ is fulfilled within 
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FIG. 13: Evolution with respect to for the generalized gir shift at a quark separation of |bT| — 3a — 0.36 fm. Figure 13a 
shows the results obtained at = 518 MeV for the SIDIS and DY limit |77||w| — >■ cxd. The solid data points correspond to the 
full result and empty symbols to the result obtained with just A7 in the numerator. Figure 
ensembles listed in Table H) 
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displays the full results for all 



statistical uncertainties. In contrast, for the T-even quantities we observe little systematic dependence on the staple 
direction and extent. As the staple extents are increased, our data appear to approach plateaus. Averages of the 
plateau values then provide estimates for the limit of infinite staple extents, ri\v\ — >■ ±00, which is formally required 
for all phenomenologically relevant TMDs. 

The physical length scale beyond which the influence of the gauge link extent diminishes is of the order of 0.4 fm 
for all cases considered. This observation invites speculations as to the physical background of this scale, e.g., an 
interpretation as color correlation length. The scale might also be related to a mass gap in the spectrum. If v is 
interpreted as the Euclidean time direction, the legs of the staple-shaped gauge link resemble static quark propagators. 
Considering our three-point function in this rotated frame of reference suggests that the plateau region is reached 
when the propagation time \riv\ of the static quark pair is large enough to suppress contributions from excited states 
sufficiently. 

Our numerical extrapolations to infinite staple extents, r]\v\ — >■ ±00, represent first predictions for the signs and 
approximate sizes of the generalized transverse shifts from lattice QCD. In particular, we find strong indications that 
the T-odd Sivers and Boer-Mulders TMDs are both sizeable and negative for the isovector, u — d quark combination 
in the case of SIDIS. 

Within statistical errors, we do not observe any clear trend in the data for the transverse shifts as functions 
of the CoUins-Soper evolution parameter ^ in the range ^ ~ 0...0.8. For the T-even generalized tensor charge, 
which shows a much better signal-to-noise ratio and less scatter of the data points, we can tentatively conclude that 
it is approximately constant in ^ for the accessible parameter ranges. We stress, however, that more quantitative 
predictions with respect to phenomenological analyses of SIDIS and DY experiments on the basis of QCD factorization 
will require much larger CoUins-Soper parameters C S> 1. For the TMD ratios discussed in this study, large C can 
in principle be accessed through larger nucleon momenta. In practice, this represents a considerable challenge due 
to quickly decreasing signal-to-noise ratios and potentially significant finite volume effects at higher P. Still, we 
expect that future lattice results for an extended range of momenta and with improved statistics will be very useful 
to establish trends in eventually allowing extrapolation into the region where factorization theorems and related 
evolution equations are applicable. 
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Appendix A: Conventions and definitions 

Whenever the four-vector b fulfills < 0, we shall make use of the abbreviation |6| = \J —b'^ . 
In the continuum, a "gauge link" or "Wilson line" is given by the path-ordered exponential 



1 



= V exp l^-ig J d\ A{Cbi\)) ' CbW j ■ (Al) 

Here the path is specified by a continuous, piecewise differentiable function Cb with derivative Ch and with C(,{0) — 0, 
Cb(l) = b. Straight Wilson lines between two points x and y shall be denoted U[x,y] and concatenations of several 
straight Wilson lines (i.e., polygons) U[x, y]U[y, z] - ■ ■ shall be abbreviated U[x, y,z, . . .]. 

For an arbitrary four-vector w, we introduce light cone coordinates w'^ — {w^ + w^)/^/2, — {w'^ — w'^)/\/2 and 
the transverse projection w_\^ = (Q^w^ ^w^ which can also be represented as a Euclidean two-component vector 
iot = {wi,W2) = {w^,w'^), W'x-W'x > 0. The basis vectors corresponding to the -I- and — components shall be denoted 
n and n, respectively, and fulfill n-n = 1. The nucleon moving in z-direction has momentum P = P+n + (m^/2P+)?i 
and spin S = A(P+/mAr)n — A(TOjv/2P+)n -I- Sj_, = —1. We use the convention e^^^'^ = 1 for the totally 
antisymmetric Levi-Civita symbol, and introduce eij = e such that ei2 = 1. 

Appendix B: Symmetry transformation properties of the correlator 

The symmetry transformation properties of 3> used in Refs. [51] [55] need to be generalized to arbitrary link directions 
V to arrive at the parametrization of Ref. |24j . The transformation properties of the corresponding 6-dependent 



correlator $ with the gauge link (14) have already been discussed in Ref. [IH] and are restated here for completeness: 



^Eubt. (b, P, rjv) = $ltrb'^ ' (Afc, AP, A^, rjAv) , (Bl) 
5 Eubtr. (b, P, S,vv)^ $L;Tum1.. (b, P, -S, vv) , (B2) 
^Ssubt. (b. P, S, vv)] * = ^LtlC."'"'' i-b, P. S, ~vv) , (B3) 

'^Ssubt. (b: P. s, vv)] * = ^L;TI;".' (-^ p, s, vv) . m 

In the equations above, the matrices A and A.-^^^ describe Lorentz transformations of vectors a;^ — >■ A^^x'^ and 
spinors ip — >■ A-^^^ip, respectively. For any Minkowski vector w = {w^,w), the space-inverted vector is defined as 
w = {w'^, —w) . 
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